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Abstract —This paper proposes a methodology for optimiz¬ 
ing a frequency-hopping network that uses continuous-phase 
frequency-shift keying and adaptive capacity-approaching chan¬ 
nel coding. The optimization takes into account the spatial 
distribution of the interfering mobiles, Nakagami fading, and 
lognormal shadowing. It includes the effects of both co-channel 
interference and adjacent-channel interference, which arises due 
to spectral-splatter effects. The average network performance 
depends on the choice of the modulation index, the number of 
frequency-hopping channels, and the fractional in-band power, 
which are assumed to be fixed network parameters. The per¬ 
formance of a given transmission depends on the code rate, 
which is adapted in response to the expected interference to meet 
a constraint on outage probability. The optimization proceeds 
by choosing a set of fixed network parameters, drawing the 
interferers from the spatiaf distribution, and determining the 
maximum rate that satisfies the outage constraint. The process is 
repeated for a large number of network realizations, and the fixed 
network parameters that maximize the area spectrai efficiency 
are identified. 

I. Introduction 

Frequency-hopping (FH) spread spectrum is a versatile 
channel access strategy for ad hoc networks m. With FH, 
a transmitter and its receiver periodically hop within a large 
set of frequency channels over a given bandwidth according 
to a pseudo-random sequence. When the transmitter and the 
receiver remain tuned to the same frequency, they are able 
to communicate. FH can exploit a much wider bandwidth 
than is possible with a single carrier frequency, and the band 
can be divided into noncontiguous segments. Furthermore, FH 
provides inherent resistance to the near-far problem without 
requiring power control. 

Continuous-phase frequency-shift keying (CPFSK) is a 
commonly used modulation for FH systems, since it has a very 
compact spectrum and produces little spectral splatter. CPFSK 
modulation is characterized by its modulation order, i.e., the 
number of possible tones, and by its modulation index h, 
which is the normalized tone spacing. For a fixed modulation 
order, the selection of h involves a tradeoff between bandwidth 
and performance. The bandwidth efficiency generally increases 
with decreasing h, while the error rate generally decreases with 
increasing h. In an ad hoc network, the tradeoff can be quan¬ 
tified by establishing a threshold on the signal-to-interference- 
and-noise ratio (SINR), which quantifies the minimum signal 
quality required for reliable communication. When a capacity- 
approaching code is used, there is a direct connection between 
the SINR threshold and the achievable rate. Together, the 


choice of modulation order, modulation index, and code rate 
constitute the modulation and coding scheme (MCS) for the 
system. 

In ||2l, the capacity of noncoherently detected CPFSK is 
computed in the absence of interference in both additive white 
Gaussian noise (AWGN) and Rayleigh fading channels. Lever¬ 
aging these results, lO optimizes a FH system in the presence 
of partial-band jamming and multiple-access interference by 
constraining the bandwidth, selecting the modulation index 
and code rate that is optimal for a Rayleigh fading channel, 
and fixing the number of frequency channels according to the 
bandwidth constraint and MCS parameters. The optimization 
is improved in H by jointly optimizing the MCS parameters 
and the number of frequency channels when the interferers are 
drawn from an arbitrary spatial distribution and the channel is 
subject to Nakagami fading and lognormal shadowing. 

The optimization of ID leverages the analysis of fSl, which 
computes the outage probability conditioned on the network 
geometry, then averages the outage probability with respect to 
the interferers’ spatial distribution. The optimization is with 
respect to the transmission capacity ii, Q, a measure of 
area spectral efficiency, which is adapted to account for the 
modulation constraints; here, we refer to it as modulation- 
constrained area spectral efficiency (MASE). While a, 11 
only consider the influence of co-channel interference (CCI), 
adjacent-channel interference (ACI) is also of concern. In H 
the analysis and optimization is extended to account for ACI 
with the goal of determining the proper amount of spectral 
splatter, which may be quantified by the fractional in-band 
power, i.e., the amount of signal energy contained within 
the hopping band. Contrary to the traditional use of a 99% 
fractional in-band power, it is found that reducing the in-band 
power to 95% improves throughput despite increasing ACI, 
due to allowing for either more frequency channels, which 
reduces the likelihood of a collision, or supporting a lower-rate 
code, which provides more error protection when a collision 
does occur. 

In all of these previous works, a globally optimal set of 
parameters is found for the network. The parameters include 
the modulation index, number of frequency channels, frac¬ 
tional in-band power, and code rate. However, by fixing these 
parameters, the outage probability will vary from one network 
realization to the next. Thus, as interferers move, join, or leave 
the network, there is a potential for wide swings in the outage 
probability, which results in an unacceptably large delay jitter 


when a retransmission protocol is used. For this reason, it is 
desirable to adapt the system in accordance with the particular 
network realization in an effort to meet a constraint on outage 
probability. When adapting a FH system, it is not desirable 
to change the number of frequency-hopping channels, the 
modulation index, or the fractional in-band power, as these 
parameters should be common to the system. Flowever, the rate 
can and should be adjusted in order to meet the performance 
objectives. 

In this paper, we consider the impact of adapting the rate 
on an ad hoc FH network. While assuming that the rate is 
adapted places this work closer to what is done in practice, it 
has a significant impact on the analysis because special care 
must be taken in how the averaging is done over the spatial 
and fading distributions. When the rate is fixed, the outage 
probability can be found by averaging over both the spatial 
and fading distributions. However, when the rate is adapted, 
the outage probability is typically constrained to not exceed 
some threshold, and thus the rate must be found that meets the 
outage constraint for each network realization. In addition to 
averaging the rate over space, the distribution of rates can be 
determined and used as a measure of throughput and fairness. 

The remainder of the paper is organized as follows. Section 
Hn describes the system model, including a discussion of the 
SINK and conditional outage probability. Section Hill discusses 
rate adaptation. Section |IV] defines MASE, describes the 
procedure used for determining the combination of parameters 
that maximize the MASE, and provides optimization results 
for several channel and network models. Einally, the paper 
concludes in Section [V] 

H. Network Model 

Consider a finite ad hoc network composed of M -f 2 
mobiles: a reference receiver, a reference source transmitter 
Xq, and M interfering mobiles Xi,X m■ The coordinate 
system is chosen such that the reference receiver is at the 
origin. The interfering mobiles are constrained to lie in a finite 
region of area A. The variable Xi represents both the 
mobile and its location, and ||Xi|| is the distance from Xi 
to the receiver. 

A. SINK 

When Xi transmits a signal with an omnidirectional an¬ 
tenna, its power at the reference receiver is 

p, = ( 1 ) 

where Pi is the average received power from mobile Xi at 
a reference distance do when both fading and shadowing are 
absent, gi is the power gain due to fading, ^i is a shadowing 
coefficient, and /(d) is a path-loss function. The power gain 
due to fading is modeled as gi = a^, where Ui has a Nakagami 
distribution with parameter rrii, and gi has unit mean. In 
Rayleigh fading, rrii = 1 and gi is exponential. While the 
{gi] are independent, they are not necessarily identically 
distributed. Eurthermore, the channel from each transmitting 
mobile to the reference receiver can have a distinct Nakagami 


parameter m^. In the presence of lognormal shadowing, the 
{/i} are independent and identically distributed zero-mean 
Gaussian with standard deviation CTs dB. In the absence of 
shadowing, = 0. Eor d > do, the path-loss function is 
expressed as the attenuation power law 


fid) = 



( 2 ) 


where a > 2 is the attenuation power-law exponent, and do 
is sufficiently large that the signals are in the far field. 

Channel access is through a frequency-hopping protocol. 
The hopping is slow, with multiple symbols per hop, which 
is a more suitable strategy for ad hoc networks than fast 
hopping m. An overall spectral band of B Hz is divided into 
L contiguous frequency channels, each of bandwidth BjL 
Hz. Let D, 1 be the duty factor of ..Aj, which is the 
probability that the mobile transmits any signal, here assumed 
to be the same for each mobile {Di = D,\/i). The system 
is characterized by two types of interference: a) CCI, which 
is when the source and the interfering mobile use the same 
frequency and occurs with probability p^, and b) ACI, which 
is when the source and the interfering mobile select adjacent 
channels and occurs with probability pa- 

Each Xi independently selects its frequency channel with 
probability 1/L and transmits with probability D. Therefore, 
the probability of CCI is pc = D/L. The frequency channel 
at the edge of the band is selected by the source Xo with 
probability 2/L, while an interior channel is chosen with 
probability {L — 2)/L. A mobile Xi, i > 0, transmits by using 
a frequency channel adjacent to the one selected by the source 
Xq with probability 1 jL if Xo selected a frequency channel at 
the edge of the band, where there is only one adjacent channel; 
otherwise, the mobile selects an adjacent channel of the source 
with probability 2/T. It follows that the probability that Xi 
produces ACI is 



2D{L-l) 

L2 


(3) 


When Xi transmits, only a fraction of the total transmit¬ 
ted power lies within the occupied frequency channel. This 
fraction is called fractional in-band power, and is denoted 
by -0. Typically, 0.95 < 4> < 0.99. The remaining fraction 
of the power is assumed to spill significantly only into the 
frequency channels that are adjacent to the one selected, which 
is a reasonable assumption for most practical systems. Let 
Kg = (1 —'0)/2 be the adjacent-channel splatter ratio, which 
indicates the fraction of power that spills into each adjacent 
channel. 

By using a physical interference model a, the instanta¬ 
neous SINR at the reference receiver is 


tpgo^o 

M 

SNR ^ ^ ^ ligAdj 

i=l 


( 4 ) 







where SNR = cIqPq/M is the signal-to-noise ratio (SNR) 
when the transmitter is at unit distance and fading and shad¬ 
owing are absent, M is the noise power, 




i = 0 
i> 1 


(5) 


is the normalized received power due to and is a 
discrete random variable that may take on three values with 
the probabilities 


with probability 
li = { with probability 

with probability pn 


( 6 ) 


where Pn = 1—Pc~Pa = 1 — DiZL — 2)/L? is the probability 
of no collision. The three probabilities correspond to the case 
that Xi selects the same frequency channel as Xq, Xi selects 
a frequency channel that is adjacent to the one chosen by 
Xq, or Xi selects a distant frequency channel (more than two 
channels from the one chosen by Xq). ACI can be neglected 
by setting ip = 1 and Ks = 0, in which case the model 
specializes to the one in i). 

B. Conditional Outage Probability 

Let /3 denote the minimum SINR required for reliable re¬ 
ception and VI = {rioi ^m} represent the set of normalized 
received powers. An outage occurs when the SINR falls below 
/3. When conditioned on VI, the outage probability is 


e(ll) = P [7 < 


(7) 


Assuming that the Nakagami parameter mo of the source’s 
channel is integer-valued and under the assumption of block 
fading, the analysis of Q using (|6l) gives 


e{Vl) = exp ( — 




mo —1 


/ /3o 


SNR/ ^ VSNR 
where /3o = /3mo/ and 


Htm 


t=o 


sm\s-t)\ 


M 


HtiVl) = 


E 
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(9) 


i>0 i=l 




The summation in (|9]l is over all sets of positive indices that 
sum to t, and from lH] 

Gui^i) =Pn5u + N +PB(pi{Ks)\ (10) 

r(L -f l)r(mj) 

where 5t, is the Kronecker delta function, equal to 1 when 
f = 0, and zero otherwise, T (cc) is the Gamma function, and 

^ + a.) 

\mi J \ rrii ) 

The conditional outage probability is calculated once the lo¬ 
cation of the mobiles and the shadowing factors are specified. 
Thus, performance can be calculated for any snapshot of the 
network topology. 


III. Rate Adaptation 

In FH systems, the probability of a collision can be reduced 
by increasing the number of frequency channels, and therefore 
it is highly desirable to choose a data modulation that has 
a compact spectrum and a reduced spectral splatter. Because 
CPFSK has these characteristics III , it is the preferred modula¬ 
tion for FH systems. When used for FH applications, CPFSK is 
typically detected noncoherently since rapid hopping precludes 
carrier synchronization. A CPFSK modulation is characterized 
by its modulation order q, which is the number of possible 
tones, and by its modulation index h, which is the normalized 
tone spacing. For FH networks, performance also depends on 
the number of frequency channels, the rate R of the error- 
control code, and the fractional in-band power. 

The rate is often adapted (e.g., cni, HD) by assigning 
a different code rate based on the quality of the channel 
to increase the system capacity and/or to satisfy a specific 
outage constraint. While in IS), an optimal common rate is 
found together with the modulation characteristics in order to 
maximize the throughput of a FH system, in this work the 
effect of link adaptation is taken into account by assuming 
that a mobile is able to estimate the channel and appropriately 
adapt its rate in order to meet a given outage constraint. 
The adaptation assumes that the distribution of the outage 
probability is known; i.e., that the transmitter has some way 
to estimate © for its intended receiver. This function can 
be precisely estimated if the transmitter locations, transmitter 
powers, shadowing coefficients, and the parameters in (|6]l, i.e. 
p). Kg, Pc and p^, are known. Importantly, the instantaneous 
fading gains or the particular channel selected by the transmit¬ 
ter need not be known. In the absence of these parameters, a 
simple adaptive scheme could request a higher rate when the 
observed outage probability (as determined by a CRC check) 
is below the target and requests a lower rate when it is above. 

Let < 7 ( 7 ) be the capacity as a function of the SINR 7 . 
The channel is in an outage when ( 7 ( 7 ) < R. Assuming 
Gaussian interference, the channel is conditionally Gaussian 
during a particular hop and the modulation-constrained AWGN 
capacity can be used for < 7 ( 7 ). The maximum achievable rate 
of noncoherent CPFSK is given in El for various modulation 
indices h, where it is called the symmetric information rate. 
In particular. Fig. 1 of El shows the symmetric information 
rate of binary CPFSK as a function of 7 for various h. To 
emphasize the dependency of the capacity on h, C(h, 7 ) 
is used in the sequel to denote the rate of CPFSK with 
modulation index h. For any value of h and 7 , the achievable 
code rate is /? = C{h,y) for the given 7 = /3. When link 
adaptation is used, the SINR threshold /3 is selected such that 
P\l ^ where e is a constraint on the channel outage. 

The appropriate /3 varies according to the number of frequency 
channels available and the fractional in-band power chosen. 

IV. Network Optimization 

For an ad hoc FH network that uses CPFSK modulation and 
adaptive-rate coding, there will be an optimal set of {L, h, ip) 
for the particular spatial distribution. The optimization is 







performed with respect to the MASE, which is defined in 
subsection IIV-AI The optimization algorithm is presented in 
subsection IIV-BI and the results are given in subsection IIV-CI 


A. Objective 

The transmission capacity 0 , 13 is a measure of the area 
spectral efficiency; i.e., the maximum rate of successful data 
transmission per unit area. The transmission capacity may be 
found by multiplying the average throughput of a typical link 
by the density of transmitters 


T (A) = AE [T] 


( 12 ) 


where X = M/A and E [T] is the throughput averaged over 
the network topologies. When accounting for modulation and 
coding, the maximum data transmission rate is determined 
by the bandwidth BjL of a frequency channel, the spectral 
efficiency of the modulation, and the code rate. 

Let r\ denote the spectral efficiency of the modulation, which 
is defined by the symbol rate divided by the 100 i/'%-power 
bandwidth of the modulation. Neglecting the switching effects 
of finite-duration hops, the spectral efficiency of CPFSK can 
be found by numerically integrating the normalized power- 
spectral density given by Equation (3.4-61) of 02 and then 
inverting the result. To emphasize the dependence of rj on h 
and Ip, we denote the spectral efficiency of CPFSK as rj {h, ip). 
When combined with a rate-i? code, the spectral efficiency of 
CPFSK becomes Rr] {h, ip) (information) bits per second per 
Hz. The data rate supported by the channel is Rrj {h, ip) B/L 
bits per second. The throughput must account for the duty 
factor D and only count correct transmissions. Hence, the 
throughput is 

y ^ Ril{h,iP)BD{l - e) 

L 

By substituting (fT3l l into (fT 2 l l. and dividing by bandwidth 
B, the normalized MASE is 


r'(A) 


r \¥.[R\ri{h,ip)D{l-e) 
'B ~ L 


(14) 


In contrast with (fTSl i. this form of transmission capacity 
explicitly takes into account the code rate R, as well as 
the spectral efficiency of the modulation rj {h, ip) and the 
equivalent number of frequency channels. 

In the literature, the transmission capacity is typically inter¬ 
preted to be the spatial density of transmissions that achieves 
a target average outage probability Q. In contrast, we place 
a constraint on the conditional outage probability; i.e., the 
outage probability of a particular network realization for which 
we are adapting the rate. Thus, the MASE developed in this 
section can be interpreted as the spatial density of transmission 
that achieves a target conditional outage probability, where the 
conditioning is with respect to the network topology. 


B. Optimization Algorithm 

For the optimization, the density of potentially interfering 
mobiles per unit area A and duty factor D have fixed values. 
For a given modulation index h, a given number of frequency 


channels, and a given fractional in-band power, the average 
code rate E [i?] can be found using a Monte Carlo approach 
as follows. Draw a realization of the network by placing 
M mobiles within an area A according to a specific spatial 
distribution. Compute the path loss from each mobile to each 
reference receiver, applying randomly generated shadowing 
factors. By setting the outage equal to the outage constraint 
e, invert the outage probability expression given by (O to 
determine the adequate SINR threshold (3. Next, find the rate 
Rn = C{h, (3) that corresponds to the SINR threshold for the 
topology and store its value. Repeat this process for a 
large number of networks N. Compute the average code rate 
as follows 

1 ^ 
n—1 

The optimization can be accomplished using an exhaustive 
search by performing the following steps; 

1) Pick a value of h. 

2) Pick a value of ip. 

3) Pick a value of L. 

4) Compute the average code rate E [i?], as described 
above, and find the spectral efficiency ri{h,ip) corre¬ 
sponding to the current h and ip. 

5) Use (fT4li to compute the normalized MASE. 

6 ) Return to step [3] until all L are considered. 

7) Return to step |2] until all ip are considered. 

8 ) Return to step [T] until all h are considered. 

The above procedure will find the t'(A) for each {L,h,ip) 
considered, and the optimal value of these parameters are the 
ones that maximize t'(A). By limiting L to be integer valued, 
the number of values is finite and an exhaustive search up to 
some maximum value is feasible. For the exhaustive search 
results presented in this section, h was quantized to a spacing 
of 0.01 over the range 0 < h < 1 and ip was quantized to a 
spacing of 0.005 over the range 0.90 <ip< 0.99. 

C. Optimization Results 

As an example, consider a network with interfering mobiles 
drawn from a binomial point process (BPP) within an annular 
region. The region has inner radius r^x = 0.25 and two outer 
radii are considered: r„et = 2 and rnet = 4. If not otherwise 
stated, the number of interfering mobiles is M = 50. The 
optimization is done by considering N =10,000 different 
network topologies. The path-loss exponent is a = 3, the duty 
factor is D = 1, the source is located at unit distance away 
from the receiver, the SNR is SNR = 10 dB, and when a 
shadowed scenario is examined, ag = 8 dB. Binary CPFSK is 
used, and the rate is adapted in order to satisfy a typical outage 
constraint of e = 0.1. Three fading models are considered: 
Rayleigh fading (m^ = 1 for all i), Nakagami fading (m^ = 4 
for all i), and mixed fading (mp = 4 and rrii = 1 for i > 1). 

Fig. □ illustrates the variability of the code rate through 
its CDF in the presence of mixed fading and shadowing. The 
figure shows three curves, corresponding to the case of a dense 
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Fig. 1. CDF of the code rate when the network is optimized. Three curves are 
provided: a) dense network (M = 50); b) moderate dense network (M = 25); 
c) sparse network (M = 5). 



Fig. 2. Optimal normalized MASE as function of the number of frequency 
channels L. For each L, the modulation index is chosen in order to maximize 
the normalized MASE, while -0 = 0.96 is adopted. A dense network is 
considered (M = 50 and rnet = 2). Top curves: mixed fading. Middle 
curves: Nakagami fading. Bottom cuiwes: Rayleigh fading. For each channel 
model both a shadowed and an unshadowed scenario are considered. 


network (M = 50), a moderately dense network (M = 25) 
and a sparse network (M = 5). The outer radius of the 
annular network is fixed to = 2. The network parameters 
are independently optimized for each network density. Fig. [T] 
shows that the CDF of the code rate is quite steep for all three 
scenarios, which emphasizes the fairness in the rate selection. 
However, the curves are steeper for dense networks than for 
sparse networks, since in a sparse network there are fewer 
interferers and therefore more variability in the SINR and code 
rate from one realization to the next. 



Fig. 3. Optimal normalized MASE as function of the modulation index h. 
For each h, the number of frequency channels is chosen in order to maximize 
the normalized MASE, while ijj = 0.96 is adopted. A dense network is 
considered {M = 50 and rnet = 2). Top curves: mixed fading. Middle 
curves: Nakagami fading. Bottom curves: Rayleigh fading. For each channel 
model both a shadowed and an unshadowed scenario are considered. 



Fig. 4. Optimal normalized MASE as function of the fractional in-band power 
Ip. For each ip, both the number of frequency channels and the modulation 
index are chosen in order to maximize the normalized MASE. A dense 
network is considered (M = 50 and rnet = 2). Top curves: mixed fading. 
Middle curves: Nakagami fading. Bottom curves: Rayleigh fading. For each 
channel model both a shadowed and an unshadowed scenario are considered. 


Fig. |2] |3l and |4] show the optimal normalized MASE 
for a dense network (M = 50 and rpet = 2, respectively) 
as function of the number of frequency channels L, the 
modulation index h, and the fractional in-band power ip. 
These three figures explore the relative importance of each 
of the three parameters of the optimization. For Figs. |2] and [3 
Ip = 0.96 is used, and for each value of the given parameter 
(L or h), the other one is chosen in order to maximize the 
normalized MASE. Eor Eig. |4]the two parameters are jointly 
optimized for each value of ip. 





































TABLE I 

Results of the Optimization for M = 50 interferers. The 

NORMALIZED MASE r' IS IN UNITS OF BPS/KHZ-m^. CHANNEL 
ABBREVIATIONS: R INDICATES RAYLEIGH FADING, N INDICATES 
NAKAGAMI fading (mi = 4, FOR ALL i), M INDICATES MIXED FADING 
(mo = 4 AND mi = 1 FOR i> 1),U INDICATES UNSHADOWED 
ENVIRONMENT, AND S INDICATES SHADOWING ((Ts = 8 DB). 


f net 

Channel 

L 

E[R] 

h 


'or>t 

sub 

2 

R/U 

315 

0.07 

0.80 

0.96 

0.79 

0.50 

N7U 

280 

0.36 

0.80 

0.96 

4.42 

2.99 

mTu 

2T9 

0.41 

0.80 

0.96 

5.00 

3.56 

R/S 

320 

0.07 

0.80 

0.96 

0.76 

0.40 

NTS 

290 

0.42 

0.80 

0.96 

4.24 

2.27 

M/S 

290 

0.38 

0.80 

0.96 

4.70 

3.02 

4 

R/U 

73 

0.05 

0.84 

0.95 

0.63 

0.32 

nTu 

80 

0.35 

0.84 

0.95 

3.74 

1.87 

mTu 

75 

0.36 

0.84 

0.95 

4.09 

1.91 

RTS 

95 

0.05 

0.84 

0.95 

0.49 

0.28 

NTS 

130 

0.40 

0.84 

0.95 

2.65 

1.75 

MTS 

100 

0.35 

0.84 

0.95 

3.03 

1.80 


The results of the optimization are shown in Table U 
For each set of parameters, the {L, h, xp) that maximize the 
normalized MASE are listed, along with the corresponding 
normalized MASE . In addition to showing the normalized 
MASE when using the optimal parameters, the normalized 
MASE is shown for a typical choice of parameters: 
iL,h,ip) = (200,0.5,0.99). 

The results shown in Table U highlight the importance of 
parameter optimization, and point out that by selecting the 
optimal parameters it is possible to improve the normalized 
MASE by a factor greater than two relative to arbitrary pa¬ 
rameter selection. Optimization results are given for different 
fading channel models and for both a shadowed and an 
unshadowed scenario. Performance is very poor for the pes¬ 
simistic assumption of Rayleigh fading, while it improves for a 
more realistic mixed fading, which implies that the source and 
destination are in the line-of-sight (LOS), while the other links 
are not. Shadowing is always detrimental, and even though it 
leads to an higher code rate, it requires a larger number of 
frequency channels. While the optimal modulation index and 
fractional in-band power vary with the density of potentially 
interfering mobiles and the distance of the transmitter from the 
receiver, similarly to HI they remain constant under different 
channel models and shadowing conditions. An increment in 
the density of potentially interfering mobiles per unit area 
leads to an increment of the normalized MASE, as well as to a 
higher code rate, number of frequency channels, and fractional 
in-band power and a lower modulation index. 


V. Conclusion 

When an EH mobile ad hoc network uses noncoherent 
CPESK and adapts its rate in order to satisfy an outage 
constraint, its performance depends on the modulation index 
h, the number of frequency channels L, and the fractional 
in-band power tp. A tradeoff among these parameters exists, 
but in practice, their values are often chosen arbitrarily. In 
this paper, the spatial distribution of the interfering mobiles 
is modeled using a point process and a physical interference 
model is adopted. Using a closed-form expression for the 
outage probability conditioned on the network geometry, the 
optimal code rate can be found, and by averaging over the 
spatial distribution, the optimal values of h, L, and ip can 
be found. Using the methodology proposed here, it is shown 
that it is possible to improve the performance in terms of area 
spectral efficiency by a factor greater than two relative to an 
arbitrary parameter selection. 
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